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ABSTRACT 

Studies of the knotty jets in both quasars and microquasars frequently make use of the ratio of the 
intensities of corresponding knots on opposite sides of the nucleus in order to infer the product of the 
intrinsic jet speed (/3jct) and the cosine of the angle the jet-axis makes with the line-of-sight (cos0), via 

the formalism /a//r — ((1 + P]ct cos6')/(l — /3jet cos 6'))^^", where a relates the intensity as a function 
of frequency as cx In the cases where /3jotcos0 is determined independently, it is found that 

the intensity ratio of a given pair of jet to counter-jet knots is over-predicted by the above formalism 
compared with the intensity ratio actually measured from radio images. As an example in the case of the 
microquasar Cygnus X-3 the original formalism predicts an intensity ratio of ^ 185, whereas the observed 
intensity ratio at one single epoch is ^ 3. Mirabel & Rodriguez (1999) have presented a refined approach 
to the original formalism which involves measuring the intensity ratio of knots when they are at equal 
angular separations from the nucleus. This method is however only applicable where there is sufficient 
time-sampling (with sufficient physical resolution) of the fading of the jet-knots so that interpolation of 
their intensities at equal distances from the nucleus is possible. This method can therefore be difficult 
to apply to microquasars and is impossible to apply to quasars. We demonstrate that inclusion of two 
indisputable physical effects: (i) the light-travel time between the knots and (ii) the simple evolution of 
the knots themselves (e.g. via adiabatic expansion) reconciles this over-prediction (in the case of Cygnus 
X-3 quoted above, an intensity ratio of ~ 3 is predicted) and renders the original formalism obsolete. 

Subject headings: methods: analytical — radiation mechanisms: non-thermal — relativity — ISM: jets 
and outflows 



1. INTRODUCTION 

Relativistic jets are observed in both quasars and micro- 
quasars, and are often seen to consist of a series of discrete 
knots moving outwards from a central nucleus, believed 
to correspond to the compact object powering the out- 
flow. Measurements of the proper motions of these knots 
are often used to constrain properties such as jet speeds 
and inclination angles, and source distance (e.g. Mirabel 
& Rodriguez 1999; Hjellming & Rupen 1995). The ratio of 
the intensities of approaching and receding knots (if there 
is sufficient spatial resolution that these can be accurately 
identified) have been used (e.g. Saripalli et al. 1997) to 
constrain their Lorentz factors, via 
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where f3 = v/c is the jet speed, 9 is the inchnation angle 
of the jet axis to the line of sight, a is the spectral index 
of the emission (defined by Si, oc where Si, is the flux 
density at frequency ly) , S'app and S'roc are the flux densities 
of a corresponding pair of approaching and receding knots, 
/Xapp and ^roc are their proper motions, and fc = 3 for a 
jet composed of discrete ejecta. 

The luminosities L{t) of the knots change with time t, 
as the knots expand and the magnetic field, and hence the 
synchrotron emissivity, decreases. Thus the true flux ratio 
is 
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where tapp and t^.^c are the times at which light leaves 
the approaching and receding knots respectively in order 
to arrive at the telescope at the same time. Unless the jet 
axis is perpendicular to the line of sight however, the light- 
travel time between approaching and receding knots will 
mean that we see the receding jet as it was at an earlier 
time, when it was more compact and hence intrinsically 
brighter (but also dimmed in the observer's frame by its 
recessional motion, taken into account by the original for- 
malism), compared with the approaching jet seen at the 
same telescope time. To account for this effect, Mirabel 
& Rodriguez (1999) proposed that the flux densities used 
to calculate the ratio should be measured at equal angu- 
lar separations from the nucleus. This cannot always be 
implemented in practice however, since this will require 
interpolation unless good temporal coverage of the jets is 
available, or unless the jet is a continuous flow, in which 
case the motion of individual knots cannot be tracked in 
any case. At early times it may also be difficult to sepa- 
rate the emission from moving jet knots and a fading core 
if there is insufficient spatial resolution. Moreover, as a 
result of opacity or the presence of a broken power law, 
interpolation of the spectrum may not be straightforward 
if in the observer's frame we sample different parts of the 
spectrum at any given frequency. 

In this Letter, we present a method of using the flux 
ratios from a single image of a source to constrain the jet 
speeds without resorting to interpolation via the Mirabel 
& Rodriguez method. 
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2. FLUX RATIOS 

2.1. Simple Scalings 

A synchrotron-emitting plasnion where the particles un- 
dergo adiabatic expansion will have a power law decay in 
intensity, L{t) oc , in which case equation 2 becomes 
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and this scaling will apply to any process which gives a 
power law decay in intensity. We consider symmetric ap- 
proaching and receding jets, in which case after ejection 
at t = 0, the epochs at which photons leave correspond- 
ing points of the front and back plasmons, t^pp and tree 
respectively, are related by 
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So in this simple case equation 3 becomes 
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While this ratio is applicable to any process that gives 
a power law decay, in the adiabatically expanding syn- 
chrotron case the parameters a and ( are not independent 
so that the determination of the flux ratios by equation 5 
is actually not introducing an extra parameter. 

2.2. Optically thin synchrotron emission and adiabatic 

expansion 

The total synchrotron cmissivity from a single, optically 
thin jet knot scales as (e.g. Longair 1994) 



J{u) oc B3/2iV (7)72^.-1/2, 



(6) 



where B is the magnetic field strength, 7 is the Lorentz 
factor of an individual electron assumed to be radiating at 
a single frequency 
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and N{pi) is the total number of electrons with energies in 
the range (7, 7-1-^7) in the plasmon, given by 

7V(7,to)=A7-f, (8) 

where p is the electron index, t{) is some arbitrary reference 
time, and A is the normalisation constant. As the plasmon 
expands, nonrelativistically, from a radius Ro{to) to R{t) 
the electron energy scales as 



7 = -^70, 
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if synchrotron losses are negligible. The spectrum then 
evolves according to 
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Putting all of the above together we find that the syn- 
chrotron emissivity of an expanding plasmon is given by 

J{V) oc j,(1-p)/2b(1+p)/2^1-P. (u) 

As the plasmon expands the magnetic field strength will 
decrease and, in the case of a tangled field, we have B oc 



1, so the plasmon emissivity has a simple dependence 
on frequency and size given by 

J{V) oc ^(l-P)/2ij(l-3p)/2_ (12) 

The ratio of flux densities as seen by the observer is then 
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Although we could take the expansion of the plasmon to be 

of the form _R oc , it is particularly instructive to look at 
the case of linear expansion, 77 = 1, for which equation 13 
becomes 



-"app 



fS'r, 



1 + (3 cose 
1 -/3cos6i 



k—p 



(14) 



This is the flux ratio observed at a given instant by the 
telescope as opposed to the interpolated fliDc at equal an- 
gular separations. As a simple generic case, emission from 
a jet composed of discrete ejecta (fc = 3) from a spectrum 
of electron index p = 2 will give an exponent of unity 
for equation 14. By way of contrast, obtaining an inter- 
polated estimate at equal angular separations will, in this 
case, give an exponent of A; -h a = 3.5. The flux ratios to 
be measured in the two cases would, however, differ, being 
measured in a single image in the former case and at equal 
angular separations in the latter. A comparison of both 
methods in any given source would of course be a useful 
means of inferring a possible asymmetry in the approach- 
ing and receding jets (Atoyan & Aharonian 1997). 

2.3. Synchrotron Self- Absorption and Spectral Breaks 

When the particle spectrum contains a break or a turnover 
we must adapt the above discussion. For example, in 
Cygnus X-3 (Miller- Jones et al. (2004)) we observe two 
discrete knots, one on each side of the central nucleus, with 
evidence for a turnover due to synchrotron self-absorption. 
In this case the spectrum will have the form oc v^^^ be- 
low the turnover frequency i^q when the knot radius is Rq 
while above this frequency the spectrum is optically thin, 
Jj/ oc i/-(p-i)/2. As the knot expands to a radius R its 
emission will then take the self-absorbed form 
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while in the optically thin regime the intensity scales like 
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The critical frequency beyond which the emissivity be- 
comes optically thin is determined by 
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and the emission at that frequency, which is the peak of 
the knot spectrum, becomes 

-5p/{p+4) 

(18) 
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Turning now to the flux ratios observed at a given in- 
stant and frequency it is clear that at suflaciently early and 
late times we will have two extremes. In the former case. 
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when the observed emission from each knot is optically 
thick we will see a oc I&l'^v^l'^ spectrum from each 
and the flux ratio exponent is A; = 3 for discrete ejecta. 
However, it may be difficult to observe actual knots in 
this regime without mixing in possible nuclear emission. 
At observed frequency v the emission will become opti- 
cally thin from the approaching knot when its radius is R\ 
which is determined by 

//? \ -(3P+4)/(p+4) 

v = (\^^co^B)vA-^\ , (19) 

while the emission from the receding knot will remain op- 
tically thick, at this frequency, until the approaching knot 
has a radius of R2 which can be easily shown to be 

\ 1 — p cos J 

Therefore, when the approaching knot has a radius R 
satisfying Ri < R < R2 the observed emission at fre- 
quency v will be a mix of Doppler boosted, optically thin 
emission from the forward knot and optically thick emis- 
sion from the receding component. The flux ratio in this 
regime is now dependent on time, i.e. knot radius, and is 
given by 
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At R = R2 all of the emission becomes optically thin at 
this frequency and the flux ratios predicted by equations 
14 and 21 are equal. Therefore, the flux ratio exponent 
drops from a value of fc to fc — p as the front and then 
the receding knot become optically thin. During this time 
the spectrum at frequency u also evolves from i/^/^ to 
and the forward knot expands by a factor R2/R1, where 
both radii are frequency dependent. The time taken for 
this expansion is determined by the expansion velocity V 
of the knot. 

In reality however, it is unlikely that the turnover in 
the spectrum would occur at one single frequency. There 
would be a finite turnover region in which the spectrum 
evolved from a z/^/^ power-law to Depending on 

the width of the turnover region and the value of f3 cos 9, 
the receding knot could be in the turnover region of the 
spectrum by the time the approaching knot had become 
optically thin. In this case, the above results would not be 
strictly applicable. 

Nonetheless, these general points apply to any source 
of opacity which changes the spectral shape or indeed to 
any broken power law that might be attributable to the 
acceleration mechanism. It presents the possibility that 
the evolution of the spectrum from flares in microquasars 
may well be influenced by the light-travel time differences 
between approaching and receding knots, as outlined in 
MiUer- Jones et al. (2004). 

2.4. Caveats 

Care should be taken if the expansion mode of the plas- 

mons changes prior to the observation from which the flux 
ratio is derived. In such a case, for example a transi- 
tion from slowed to free expansion (Hjellming & Johnston 
1988), the time decay of the flux density would change 
(steepen with time in this case). In order to use flux ra- 
tios to constrain the value of (3 cos 6, the flux densities of 



the approaching and receding knots would then have to 
be measured when the knots were both in the same ex- 
pansion regime. Unless the transition radius were known, 
this would require actually measuring (as opposed to in- 
terpolating) the flux densities at equal angular separation 
from the core. We note that if there is significant deceler- 
ation of the expanding plasmons due to interaction with 
surrounding material, as mentioned by Hjellming & Han 
(1995), then (R/Rq) oc (t/to)'^ , where rj < 1, and equation 

14 then requires modification. We also draw attention to 
Fender (2003), which presents caveats to be considered 
when using proper motions to place limits on the bulk 
Lorentz factors of jets; any Lorentz factors thus derived 
are strictly only lower limits. 

3. COMPARISON WITH OBSERVATIONS 

The VLBA observations of Cygnus X-3 presented by 
Miller-Jones et al. (2004) show a jet which at 5 GHz and 

15 GHz appears to be composed of two separating discrete 
knots, which were interpreted as approaching and reced- 
ing plasmons. A precession modelling analysis yielded a 
value (3 cos 6 = 0.62 ±0.11, and the spectral index of the 
emission was found to be a = 0.60 ±0.05. Assuming linear 
expansion of the jet knots, we would thus predict a flux 
density ratio of 3.2 ± 1.0. For the last two epochs (2001 
September 20 and 21), the measured flux ratios are given 
in Table 1 . While not matching the theoretical prediction 
perfectly, they are now of the correct order, in contrast 
with the predictions of the original formalism, which is 
wrong by two orders of magnitude. There are various pos- 
sible explanations for the slight discrepancy. Most impor- 
tantly, the measurement of the flux densities themselves 
was often difficult. It is also possible that the plasmon 
expansion was not exactly linear with time, which would 
alter the exponent in equation 14 and change the predicted 
fiux density ratio. Moreover, the measured spectral index 
a was for the integrated spectrum; the values of a and p 
could in principle differ for the individual jet knots. The 
quality of the data makes it difficult to interpolate back to 
the flux densities at equal angular separations in this case, 
but our best attempts gave flux ratios between 1.58 and 
10.63. In such cases, our direct measurement method gives 
a miich more accurate determination of the expected flux 
ratio, for more meaningful comparison with the jet speeds 
and inclination angle found by different methods. We note 
that if the spectral index of the jet material is known, 
our method requires only a single image to determine the 
value of /? cos 9, whereas the interpolation method requires 
at least two images taken at different times. This frees it 
from the uncertainty inherent in comparing VLBI images, 
particularly if the imaging is difficult, as was the case in 
these observations (Miller- Jones et al. 2004). For a sin- 
gle image, the ratio of two flux densities is set, whereas 
when comparing different images, in order to be able to 
interpolate accurately, one has to be confldent that one 
has recovered the same fraction of the true flux density in 
both images in order to be able to take an accurate flux 
density ratio. 

This theory could also be applied to the observations of 
GRS 1915+105 detailed by Mirabel & Rodriguez (1994). 
They observcxl discrete radio ejecta moving outward from 
the nucleus over a period of ~ 1 month. Again, we take 
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the flux density ratio of their observed knots once they 
had clearly separated from one another and from the nu- 
cleus, and we only compare corresponding pairs of ejecta. 
From their derived value of (3 cos 9 = 0.323 ± 0.016 and 
their quoted spectral index of a = 0.84 ± 0.03, we predict 
a flux density ratio of 1.24 ± 0.05. For the later epochs 
(1994 April 16, 23 and 30), the measured flux density ra- 
tios are 2.33, 2.63 and 1.80 respectively. Again, this is 
slightly greater than we predict, but is of the right order. 
Underpredicting the flux density ratio implies the expo- 
nent should be larger in equation 14, which requires r] > 1, 
i.e. the expansion scales slightly more rapidly with time 
than Rent. 

4. CONCLUSIONS 

We have considered the evolution of synchrotron bub- 
bles (plasmons) in oppositely-directed microquasar jets. 
We have found that our new formalism can explain the 
observed flux density ratios in microquasar jets in systems 
in which the synchrotron bubble model is applicable, such 
as Cygnus X-3. In contrast, the original formalism con- 
siderably overpredicts the observed flux density ratio in 
observations of this system. In the case of free (linear) 
expansion, (R/Rq) cx (t/ta), we foimd that the flux ratios 
of the approaching and receding plasmons are given by 

^app/^rec = ((1 + /3 COS 0)/(l - /3 COS 0))'=-?' . 
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Table 1 

Cygnus X-3 VLBA flux density ratios measured from observations of 2001 September outburst of Cygnus X-3 



Date (UT) 


Observing Frequency (GHz) 


Flux density ratio (South/North) 


September 20 


5 


1.43 ± 0.05 


September 21 


5 


2.39 ± 0.10 


September 20 


15 


1.14 ± 0.19 


September 21 


15 


3.09 ±0.14 



